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Abstract 

The time-dependent susceptibility for the finite-size mean-field Ran- 
dom Orthogonal model (ROM) is studied numerically for temperatures 
above the mode-coupling temperature. The results show that the imag- 
inary part of the susceptibility x"iy) obeys the scaling form proposed 
for glass-forming liquids with the peak frequency decreasesing as the 
temperature is lowered consistently with the Vogel-Fulcher law with a 
critical temperature remarkably close to the known critical tempera- 
ture T c of the model where the configurational entropy vanishes. 

PACS: 64.40.-i, 64.60.Cn, 75.10.Nr 

The spectral properties of the primary or a-relaxation in supercooled 
liquids has been largely studied by means of dielectric spectroscopy (Dixon 
et al. 1990, Chamberlin 1991, Dixon el al. 1991, Schonhals et al. 1991, 
Menon and Nagel 1993, Schonhals et al. 1993, Kuldlik et al. 1995, Leheny 
et al. 1996, Leheny and Nagel 1997) finding that the data for the imaginary 
part e"(u) of the dielectric susceptibility e(u) at different temperatures and 
for several glass-forming liquids can be collapsed onto a master curve using 
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a three-parameter scaling function. The master plot is able to reproduce 
the e"(u) data around the relaxation peak u p and also at higher frequencies. 
The bad collapse in the low frequency part has been object of some debates 
(Schbnhals et al. 1991, Menon and Nagel 93, Schonhals et al. 1991, Kuldlik 
et al. 1995, Leheny et al. 1996), however there is no dispute above u p . The 
frequency v v has a very strong temperature dependence, commonly fitted by 
a Vogel-Fulcher form log 10 (f p ) = log 10 (fo) — A/(T — To), where To is close 
to the Kauzmann temperature (Kauzmann 1948) where the configurational 
entropy vanishes (see e.g., Angell 1988). 

In this contribution we compare the frequency-domain analysis of the 
finite-size Random Orthogonal Model (ROM) above the mode coupling tem- 
perature with the above scenario. The main motivation for this study was 
to make a stringent test on the ROM as a possible toy-model for the fragile- 
glass scenario. Model Hamiltonians capable of describing relaxation pro- 
cesses in supercooled liquids and structural glasses are difficult to obtain. 
However, starting with the work of Kirkpatrick, Thirumalai and Wolynes 
(Kirkpatrick and Thirumalai 1987a, b, Kirkpatrick and Wolynes 1987) in 
the late 80's, it is now clear that there is a close analogy between some 
mean- field spin-glass models and structural glasses (Bouchaud et al. 1998). 
The basic simplification occurring in mean-field models is that in the limit 
of a very large (N — > oo) number of spin one is left with a closed set of 
equations for the two-time correlation and response functions which, above 
a critical temperature Tp, are equivalent to the schematic mode coupling 
equations introduced by Leutheusser, Gotze and others (Bengtzelius et al. 
1984, Leutheusser 1984, Gotze 1991) as a model for the ideal glass transition. 

In mean-field models the barrier separating different ergodic components 
diverges in the mean-field limit, hence at the critical temperature Tjj a real 
ergodic to non-ergodic transition takes place with diverging relaxation times. 
The critical temperature Tp coincides with the critical temperature Tmct 
derived in the Mode-Coupling Theory (MCT) and in what follows we will 
use only the notation Tmct- In a real system, however, barriers are of 
finite height and the glass transition appears at T g < Tmct where the 
typical activation time over barriers is of the same order of the observation 
time. In these systems Tmct represents the temperature below which the 
dynamics is dominated by activated hopping over energy barries. Therefore 
to go beyond mean-field it is necessary to include activated processes, a very 
difficult task since it implies the knowledge of excitations involved in the 
dynamics. Recent studies have shown that activated processes in mean-field 
models could be included just keeping N finite (Crisanti and Ritort 2000a, 
b), giving support to the scenario of the fragile glass transition developed 
from spin-glass models. This is not a trivial assumption since it is not a 
priori clear why excitations in mean-field spin glass models should have 
similar properties to those of supercooled liquids. This, for example, seems 
to be the case for the Random Orthogonal model (ROM) (Marinari et al. 
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1994), but not for the mean- field Potts-Glass model (Bragian et al. 2001, 
Bragian et al. 2002). 

To compare the results from the finite-size ROM with the experimental 
one (Dixon et al. 1990, Chamberlin 1991, Schonhals et al. 1991, Dixon 
et al. 1991, Menon and Nagel 1993, Schonhals et al. 1993, Kuldlik et al. 
1995, Menon and Nagel 1995, Leheny et al. 1996, Leheny and Nagel 1997) in 
this contribution we shall consider only temperatures above Tmct- Since the 
range of temperatures we explore are all above the mode-coupling transition 
Tmct, we do not expect to find diverging timescales in the large N limit. 

The ROM (Marinari et al. 1994) is defined by the Hamiltonian 

H = -2^2 Jij cj,i cfj - °i (!) 

ij i 

where <7j = ±1 are A" Ising spin variables, and J™ is a N x N random 
symmetric orthogonal matrix with Ju = 0. For N — > oo and h = this 
model has a dynamical transition at Tmct = 0.536, and a static transition 
at T c = 0.256... (Marinari et al. 1994). Numerical simulations are performed 
using the Monte Carlo (MC) method with the Glauber algorithm for tem- 
peratures in the range 0.6 up to 2.0. To study the frequency response we 
considered a time-dependent field of the form h(t) = ho cos(2 ir ut), where 
the time is measured in MC steps and h,Q = 0.2 small enough to be within the 
linear response regime. In our simulations the typical range of v was 10~ 6 - 
10 _1 . For each frequency v the complex susceptibility xiy) = x'(y) +^/t"( zy ) 
is given by 

, M N 

= WmT.T,^) COB(27TI/t), (2) 
t=l j=l 

, M N 

x» = t^EE^w^ 2 ^*)- ( 3 ) 

The number of MC steps M after equilibration was 100 for the largest v 
and up to 10 7 for the shortest v. As system size we used = 300 which is 
a good compromise between small sample-to-sample fluctuations and small 
barriers height. 

Figure ^ shows the real and imaginary parts of the susceptibility over 
the available range of frequency. Not all temperatures are reported for a 
better drawing. The relaxation peak in the imaginary part can be fitted 
with a log-normal form (Wu and Nagel 1992): 

X» = exp \-(log w u - log^VS 2 ! (4) 

where u p is the frequency of the peak, E the width, and Ax = x'o ~ X'ooi 
where x'o an d x'oo are > respectively, the low and high frequency limit of x'( u )- 
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Figure 1: The real and imaginary part of the complex susceptibility x' an d 
x" as function of v for the ROM with N = 300 at different temperatures. 
Upper panel top to bottom, lower panel left to right, T = 0.7, 0.9, 1.1 and 
1.3. 

As the temperature is lowered the peak frequency u p decreases, and the 
width S broadens. The behavior of v„ is consistent with the Vogel-Fulcher 
law exp[— A/(T — To)] (Menon and Nagel 1995). The fit of the frequency 
peak Up for the ROM with the Vogel-Fulcher formula is rather good, see Fig. 
El and gives A = 0.89 ±0.06 and lnz/ = 0.64 ±0.02 T = 0.28 ±0.02, a value 
in agreement with the critical value T c = 0.256.... We note, however, that 
data can also be fitted using different expressions such as exp[— A/T 2 ] or 
the Adam-Gibbs formula exp[— A/(TS C (T))] where S c is the configurational 
entropy, and in particular with the formula (CR) v p = uq exp[— A/3 e g(T)/T] 
where (3 e ft{T) = dS c (ei s ) / dei s \ e . a=e . a ^ derived from a cooperative scenario of 
relaxation (Crisanti and Ritort 2002), see Fig. |2j The CR formula predicts 
a crossover form fragile to strong behaviors as the temperature is lowered, 
however, differences among all these expressions can be appreciated only for 
very low values of v p which are out of our measurements range. 

The analysis of the response for glass-former liquids reveals three power 
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Figure 2: log 10 f p as function of T for the ROM with N = 300. The full line is 
the Vogel-Fulcher law u p = Vq exp[— A/ (T — To)] while the dashed line is the 
formula u p = vq exp[— A(3 e s(T)/T] of Ref. (Crisanti and Ritort 2002) with 
/? e ff(T) = dS c (e\s) I ' dei B \ eiB=ei ji>) evaluated using the results of Refs. (Crisanti 
and Ritort 2000a, b). The discrepancy at high temperature is probably 
due to a poor numerical estimation of the configurational entropy, indeed a 
similar deviation is found using the Adam-Gibbs formula (not reported). 

laws for x" (Leheny and Nagel 1997): 



The discrete nature of Monte Carlo dynamics time step prevents us from 
resolving the last, nevertheless the first two regimes are clearly seen, as 
shown in Figure [3 At higher temperatures m = (3 = 1 and the relaxation 
is Debye-like with exponential decaying correlations. As the temperature 
is lowered the value of (3 decreases below 1 and decay becomes stretched- 
exponential. It is known that for glass-former liquids (3 and a are related 
by (<t + l)/{f3 + 1) = 7, where 7 is a constant (Leheny and Nagel 1997). 
Furthermore a varies linearly with temperature: a = B(T—T a ) with T a ~ To 
(Leheny and Nagel 1997), This implies that (3 = B'(T - T ) + (1 - 
Inserting into this formula the values of (3 obtained for the ROM at various 
temperatures and the value of To computed from u p we find 7 = 0.72 ± 0.02 
the same value found for real liquids (Menon and Nagel 1995, Leheny and 
Nagel 1997). 

The analyticity of x{ v ) an d linearity of absorption at asymptotically low 
frequencies implies that x"{ v ) °^ v f° r v ^ v v (Schonhals et al. 1991, Menon 
and Nagel 1993, Schonhals et al. 1993). For the ROM with N = 300 we 
find m ~ 1 for temperatures down to about T = 0.8 while below significant 




v < v p 

u > Up 
v » v „ 
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Figure 3: x"W) for the R0M witn ^ = 300 for temperatures T = 0.7 and 
1.1. The exponents are m = 1, f3 = 1. for T = 1.1 and m = .96, (3 = .71 for 
T = 0.7. Inset: x"iy) & t temperature T = 0.6 for N = 64 (filled circles) and 
300 (empty circles). The lines have slope m = 0.3 for N = 64 and m = 0.5 
for N = 300. The increase of m toward 1 as N grows is clearly seen. 

deviations with m < 1 are observed. Similar deviations have been observed 
in data from glass-forming liquids and generated some controversy (Dixon 
et al. 1990, Chamberlin 1991, Schonhals et al. 1991, Dixon et al. 1991, 
Menon and Nagel 1993, Schonhals et al. 1993, Kuldlik et al. 1995, Leheny 
et al. 1996) on the reliability of the scaling form proposed by Dixon et al. 
(Dixon et al. 1990, Chamberlin 1991, Dixon et al. 1991). Many liquids 
posses secondary relaxations which overlap the primary response broaden- 
ing the peak and leading to deviation from linearity (Kuldlik et al. 1995, 
Leheny et al. 1996). In the case of ROM these secondary relaxations are 
related to the fact that the barriers separating the low states sampled as 
the temperature is decreased toward Tmct are not well separated for not 
too large N. Indeed studies of mean-field spin-glass models for the struc- 
tural glass transition shows that in the thermodynamic limit there is no gap 
between saddles separating local minima with energy above the threshold 
energy associated with the dynamical transition (Cavagna et al. 1997). This 
is a situation more reminiscent of spin-glasses rather than glasses for which 
both experimental (Bitko et al. 1996) and numerical simulations (Bitko et 
al. 1996, Rao 2001) show a broader shape of x"( v ) near the peak. 

This scenario is supported by a finite-size scaling analysis of the ROM. 
Indeed we find that for a fixed temperature while f3 is independent on N, 
the value of m, when less than 1, increases toward 1 as N is increased. 
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In the inset of Figure El we show x"(^) at T = 0.6 for systems of size 
N = 64 and N = 300. The increase of m going from 64 to 300 spins 
is clear. Finally we address the goodness of the Nagel scaling. Dixon et 
al (Dixon et al. 1990, Chamberlin 1991, Dixon et al. 1991) have shown 
that all data for the dielectric susceptibility e{v) of different glass-forming 
liquids and temperatures can be collapsed onto a single master curve by 
plotting (l/tt))log 10 (e"(i')z/p/z/Ae) versus (l/w)(l/w + l)log 10 (^ / V p ) , where 
Ae = e'(0) — e'(oo), u p is the peak frequency and w is the half-maximum 
width of the e"{v) peak normalized to the corresponding width of the Debye 
peak: wd — 1-14 decades. 




-2 
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Figure 4: The Nagel's scaling for the ROM with N = 300 and temperatures 
T = 2.0, 1.3, 1.2, 1.1, 1.0, 0.9, 0.8, 0.7 and 0.6. 

Figure |H shows the Nagel plot for the ROM with N = 300. The data are 
from temperatures ranging from 0.6 up to 2.0. For each curve the parameters 
Ax, w = 2\^2Tj/w,i and v v have been obtained from the log-normal fit of 
x" using (j3|. We see that while the collapse for v > u p is good for all 
temperatures, for u < v p only data with m = 1 do collapse. As noted in 
Refs. (Kuldlik et al. 1995, Leheny et al. 1996) the optimization of the three 
parameters through the fitting is essential to have a good collapse of data. 

In conclusion, we have shown that the primary relaxation in the finite size 
mean-field ROM obeys the scaling form typical of glass- forming liquids. Fur- 
thermore, the frequency peak of the imaginary part of the complex suscepti- 
bility follows the Vogel-Fulcher law with critical temperature To = 0.28±0.02 
very close to the critical temperature T c = 0.256..., the Kauzmann temper- 
ature of the model. All system sizes studied (up to N = 300) lead to this 
value for To. Because we used Monte Carlo dynamics there is a maximum 
value for the frequency ~ determined by the discreteness of the elemen- 
tary time step. As a consequence we are not able to resolve the second 
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scaling behaviors of the imaginary part of the susceptibility x"(f) ~ v~ a 
for v S> v p . Nevertheless by assuming that the exponent a vanishes linearly 
at To, we obtained for the constant 7 relating the exponents (3 and a the 
same value 7 = 0.72 ± 0.02, found for real glass-forming liquids (Menon and 
Nagel 1995, Leheny and Nagel 1997). This value does not depend on the 
system size for the sizes we studied. Overall, the present results show that 
finite-size mean-field spin glasses capture the cooperative effects responsi- 
ble for the relaxational processes observed in glass forming liquids when 
approaching the mode-coupling temperature from above. The extension of 
this analysis to the region below Tmct where strong finite N effects are to 
be observed remains an interesting open problem. 

We acknowledge F. Sciortino and P.Tartaglia for a critical reading of the 
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